Magnetic anisotropy of single 3d spins on CuN surface 
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First-principles calculations of the magnetic anisotropy energy for Mn- and Fe-atoms on 
CuN/Cu(001) surface are performed making use of the torque method. The easy magnetization 
direction is found to be different for Mn and Fe atoms in accord with the experiment. It is shown 
the magnetic anisotropy has a single-ion character and mainly originates from the local magnetic mo- 
ment of Mn- and Fe-atoms. The uniaxial magnetic anisotropy constants are calculated in reasonable 
agreement with the experiment. 
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Recent scanning tunneling miscropscopy (STM) mea- 
surements of the spin-excitation energies in a mag- 
netic field [l| for individual Fe and Mn atoms on 
CuN/Cu(001)-c(2x2) substrate report large values of the 
axial and transverse magnetic anisotropy energies (MAE) 
for a single magnetic atom. The STM experiments are 
complemented by density-functional theoretical calcula- 
tions. These calculations reveal that the magnetic atoms 
become incorporated into a covalent CuN matrix, so that 
their electronic and magnetic character differs from the 
gas-phase transition metal atoms. 

These STM experiments [l| along with previously re- 
ported XMCD measurements for a single Co atom 
and small Co clusters on the Pt(lll) surface show that 
just a few atom size nanostructures can maintain a stable 
magnetic orientation at low temperature due to the large 
magnetic anisotropy energy (MAE). What makes these 
atomic-scale magnetic structures technologically relevant 
is their large MAE which provides the means of reducing 
the size of the magnetic bits above the superparamag- 
netic limit, i.e. the ratio of the MAE to the thermal 
energy fcsT. Understanding of the atomic-scale MAE in 
nanomagnets is essential in the determination of the min- 
imum feasible magnetic memory bit size, and can assist 
in further increase of the magnetic recording density. 

In the work reported here we make use of ab initio nu- 
merical calculations of the MAE to analyze the key physi- 
cal quantities determining the anisotropic magnetic char- 
acteristics of single 3d-metal atoms on CuN/Cu(001)- 
c(2x2) substrate. Similar to the theory of Ref. we 
use a supercell model. The supercell consists of three 
Cu(OOl) layers and a single CU2N atomic layer with 
c(2x2)N-Cu(001) arrangement given in The in-plane 
c(2x2) dimentional unit cell is doubled (CU4N2), and the 
3d-atom (Mn and Fe) is placed on the top of Cu-atom 
. The supercell is shown schematically in Fig. 1. The 
vacuum is modeled by the equivalent of four empty Cu 
layers. 

The structure relaxation is performed employing the 
standard VASP method Q without spin-orbit coupling 
(SOC) and making use of the generalized gradient ap- 
proximation. Placing 3d atom on the top of the Cu atom 
in the CuN surface makes a substantial rearrangement of 



the atomic structure (see Fig. 1). The Cu atom right 
below the adatom moves toward the bulk and the re- 
laxed distance between this atom and magnetic atom 
is decreasing from 4.42 Bohr for the Mn atom to 4.27 
Bohr for the Fe atom. Other Cu atoms in the CuN top- 
layer change slightly their positions with the change of 
the magnetic atom. Overall relaxed atomic positions are 
qualitatively consistent with the picture given in Ref. [l| 




FIG. 1: A schematic crystal structure used to represent the 
3d-atom on the c(2x2)N-Cu(001)surface. The actual atomic 
positions correspond to the case of the Mn atom on the 
c(2x2)N-Cu(001)surface. 

We investigate the relativistic electronic and magnetic 
character of 3d-atoms on the c(2x2)N-Cu(001)surface. 
We use the relativistic version of the full-potential lin- 
earized augmented plane-wave method (FP-LAPW) H, 
in which spin-orbit (SO) coupling is included in a self- 
consistent second- variational procedure i 6] . The conven- 
tional (von Barth-Hedin) local spin-density approxima- 
tion is adopted in the calculations, which is expected to 
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FIG. 2: Spin-resolved PDOS for Mn and Fe adatoms. Also 
shown PDOS for the Fe atom — and iz^ — r'^ spin- 
minority orbitals. 

be valid for itinerant metallic systems. 

The spin Ms and orbital Ml magnetic moments for 
the magnetization directed along the z-axis are given in 
Table I. for the Mn and Fe atoms. Small spin and orbital 
moments are also induced on neighboring Cu sites and 
quickly decay away from the magnetic Mn or Fe atom. 
The spin-resolved projected density of states (PDOS) for 
the Mn and Fe atoms is shown in Fig. 2. The spin- 
majority manifold is practically fully occupied for both 
Mn and Fe. For Mn atom, the spin-minority channel 
is almost empty and the orbital Ml moment is almost 
zero. The spin-minority occupation is increased for the 
Fe-adatom while the spin splitting and spin moment Ms 
are decreasing. The detailed inspection of m;-projected 
PDOS shows that non-zero orbital M^ moment for the Fe 
atom originates from \ms = —^;mi = ±2) orbitals near 
Fermi edge. The major contribution to Ml is brought 
about mainly by in-plane xy and — spin-minority 
orbitals. The — r^(~ jm; = 0}) spin-minority orbital 
(see Fig. 2) does not contribute to Ml- This out-of-plane 
— orbital is the least localized due to the strong 
overlap with 3d electrons of the Cu atom beneath. 

Next we turn to a salient aspect of our investigation, 
the MAE calculations. The anisotropic energy EA{d,(f>) 
dependence (including the second order terms) on the 
magnetization direction reads, 

EA{e,c^) = K^el+K^el-el), (1) 
EAie, (f) = coa^{e) + k}} sm^{e){coa^{(j)) - ain^ {(/>)) 



TABLE I: Total spin moment per unit ceU (A'/s^°'), spin (Ms) 
and orbital (Ml) magnetic moments on 3d-adatom (in Bohr 
magnetons) for the magnetization directed along the z-axis. 

Atom Ms^°* Ms Ml 
Mn 4.379 3.758 0.004 
Fe 3.654 2.917 0.076 



where and K2 are the uniaxial MAE constants, and 
Sx,y,z are the cartesian coordinates of the normalized 
magnetization vector M/|]VI|. The 9 and <j> arc the polar 
angles in the reference frame which is chosen as follows: 
the X-axis is along the in-plane hollow direction, the y- 
axis is along the in-plane N-chain direction, and z-axis is 
along the out-of-plane direction (see Fig. 1). 

In order to evaluate the MAE from Eq.([T]), we make 
use of the torque method ■ It can be formulated as 
follows. We solve the Kohn-Sham equations for a two- 



component spinor 1$^) = 




J2(-y' + V,.ff+ai-s)) <i>f(r) = e.$?(r), (2) 

where the Veff = V{r)I + cr ■ B(r) matrix consists of 
the sum of the scalar potential V and "exchange" field 
B parallel to the spin moment Ms, and ^^gQ = C(l ' ^) 
is the SO coupling operator. When the magnetic force 
theorem Q is used to evaluate the magnetocrystalline 
anisotropy energy, the Ms is rotated and a single en- 
ergy band calculation is performed for the new orienta- 
tion of Ms- The MAE results from SO coupling induced 
changes in the band eigenvalues Ea{0, (p) — 4>)- 
Alternatively, the torque T{9, (p) — dEA{d, 4')ldO can be 
evaluated making use of the linear response theory: 

occ T /~)TTt 

m<^) = 5](ci>:i— c(r-s)ut+ue(r.s)— 1$:) (s) 

i 

where the U(0, (p) is a conventional spin rotation ma- 
trix and 1$') = U(6', (/))|<I>). An advantage of this ap- 
proach is that it allows the total MAE separation into 
the element-specific contributions from different atoms 
in the unit cell. The torque method has been first im- 
plemented in FP-LAPW basis in Ref. Also, it has 
been employed recently in the Korringa-Kohn-Rostocker 
calculations 

The torque T(9, (p) angular dependence is shown in 
Fig. 3 for both Mn and Fe-atoms on CuN/Cu(001). A 
set of 784 /c-points in the full 2D-BZ which is equivalent 
to 3136 /c-points in the full 2D-BZ of Cu(OOl)) is used 
in these calculations. The uniaxial MAE constants K2 

and K2 can be evaluated from the torque T(0, (p) angu- 
lar dependence, shown in Fig. 3 and angular derivative 
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dence shown in Fig. 3, we obtain MAE= J^^^ dO T{e, (j)). 
The values of the MAE are shown in Table II. There is 
an increase of the MAE from Mn to the Fe atom case. 
For the Mn atom, the easy magnetization axis is directed 
along the surface normal z-axis, in agreement with the 
experimental data For the case of Fe, the easy mag- 
netization is along the N-chain, also in agreement with 
the experiment m. The anisotropic energy EA{0,(j)) an- 
gular dependence for Mn and Fe atoms on CuN surface 
together with the easy magnetization axis orientation is 
illustrated in Fig. 4. 



Total MAE Mn Fe 

AEa[z~x] -0.03 -0.77 

AEaIz-v] -0.37 1.08 

AEaIv-x] 0.34 -1.86 

TABLE II: The MAE (meV) for the Mn and Fe-adatoms, 
Here, AE[i ^ j] = EA[M\\j] - EA[M\\i], i{j) ^x,y,z 



o 




• Fe(Phi=0) 
■ Total (Phi=0) 

* Fe (Phi=90°) 

♦ Total (Phl=90°) 



♦ 



60 80 100 120 
2*Theta (degrees) 




FIG. 3: The torque T{e, (/>) for the x - z-plane (<^ = 0°) and 
y — 2;-plane {(j) = 90°) as a function of 0. The total torque and 
leading contributions from the 3d-adatom are shown. 



of Eq.©, 



Kij cos(2(?!))]sin(26'). 



(4) 



stants are: K2- 



For the Mn atom, the values of the uniaxial MAE con- 
0.20 meV and fs:i'=-0.17 meV. The Mn 
atom contribution in K2 = -0.16 meV, and K2—-Q.\2 

meV. For the Fe atom, K^={).1& meV and K^^=Qm 
meV, and the Fe atom specific contribitions in K2 = - 

0.04 meV, and =0-97 meV. Also, we found that higher 
order anisotropy is much less (at least by an order of mag- 
nitude) than the uniaxial anisotropy. 

Now we evaluate the MAE defined as the energy 
Ea{(^,4') difference for different directions of the mag- 
netization M. Using the torque T(0, 0) angular depen- 



It is quite common to examine the correlation between 
the MAE and the orbital moment anisotropy (OMA). 
Approximate relation between the MAE and OMA is 
given by Bruno formula [MAE « -^/4 OMA], where 
^ is the SOC constant. For the Mn atom case, the Bruno 
formula gives AEa[z — x] of -0.07 meV, AEa[z — y] of 
-0.02 meV, and AEaIu — x] of -0.05 meV. Comparison 
with the torgue results of Table II. shows that Bruno for- 
mula yields the correct easy z-axis but fails to describe 
the y — X plane transverse anisotropy. For the case of Fe 
atom, making use of Bruno formula we obtain the MAE 
of AEa[z - x] = -0.74 mcV, AEa[z - y] ^ 1.46 meV, 
and AEA[y — x] = -2.19 meV in a good agreement with 
the torgue results (see Table II.). The reason why Bruno 
formula works better for Fe atom than for Mn atom case 
is that it is not accurate enough to account for relatively 
small Mn atom MAE. For the stronger Fe atom MAE, 
validity of Bruno formula is improving on qualitative. 

Now we turn to comparison with the experimental re- 
sults of Hirjibehedin et al. [1]. The STM measures the 
spin-excitation energies in a magnetic field. These excita- 
tion spectra are then analysed by the model Hamiltonian, 



H = .g/isBS + DSl + E{Sl - Sl) 



(5) 



In the Eq.(l5]), the z-axis is chosen along the easy magneti- 
zation direction. In order to compare our results with the 
experiment, we have to convert the data in Table II. into 
the reference frame chosen in Ref. [l^ and re-normalize 
the anisotropy values by 5^ {S = 5/2 for Mn and S = 2 
for Fe). The results are shown in Table III. Our ab initio 
results correctly reproduce the sign and order of magni- 
tude of D and E experimental anisotropics. 

It is quite surprising that the LSDA based calculations 
give quite reasonable values of the MAE constants for 
the systems which have been initially thought as being 
close to the atomic limit. The density functional theory 
is known to work reasonably well for the ground state 
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FIG. 4: The anisotropic energy {9, (j)} angular dependence 
for Mn (top) and Fe (bottom) atoms on CuN surface. Here 
the length of radius vector forming the surface is equal to 
[EA{d,4>) — EA{esLsy axis)]. 



energy determination, and the MAE is defined as the 



ground state energy difference for different magnetization 
directions. Most probably, that is why the MAE results 
of density functional theory resemble the values of the 
uniaxial MAE constants experimentally determined from 
the spin excitation spectra of the atomic spin. On the 
other hand, the electron correlation effects beyond those 
which are already included in LSDA can play essential 
role in more accurate theoretical modeling and interpre- 
tation of the experimental data [J| . Further progress in 
realistic calculations of the ground state properties and 
excitations for single atomic spin on the surface will be 
made on the basis of the newly emerging combination of 
the LSDA and dynamical mean field theory p^ . 

TABLE III: Comparison with experimental D and E (meV). 
Here, for the Mn atom case, the x-axis is along N-chain, the 
y-axis is along the hollow direction, and the 2-axis is out-of- 
plane. For the Fe atom case, the x-axis is along the hollow 
direction , the j/-axis is out-of-plane, and the z-axis along N- 
chain 0]. 



Mn 


D 


E 


(meV) 


Exp. 


-0.039 ± 0.001 


0.007 ± 0.001 




LSDA 


-0.03 


0.03 




Fe 


D 


E 


(meV) 


Exp. 


-L55 ± 0.01 


0.31 ± 0.01 




LSDA 


-0.36 


0.10 





In conclusion, we have shown that the mag- 
netic anisotropy energies for Mn- and Fe-atoms on 
CuN/Cu(001) can be semi-qualitatively reproduced by 
the first-principles LSDA FP-LAPW calculations. The 
easy magnetization direction is found in agreement with 
the experimental data for Mn and Fe atoms. It is shown 
the calculated MAE has a single-ion character and mainly 
originates from the well localized moment of Mn- and 
Fe-atoms. The uniaxial MAE constants are calculated in 
semi-quantitative agreement with the experiment. 
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